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1.   Introdqction. 

Kopal  [3^.  bAS  sucgttstttd  a  nfrteal  mmttned  tor   Bolrlng  dlff«r«atljd 
•quAtlona  uhloh  ums  rational  Pad*  ooirr«rc«t«  to  approxlMit*  fwotlOBa 
of  tho  forward  ihlft  oporator  I  to  caaorato  rMvrsioa  focwtl**.  Ba  baa 
daaonatratad  tha  aarlta  of  nalng  aooh  f onnlaa  la  aelTlng  flrat  Mdar 
aquatlona  and  haa  «dilblt«d  oartaln  atabilltgr  problaaa  ralatad  tiMrato* 

Iha  purpoaa  of  thla  papar  ia  to  axtand  tha  atu^  of  thla  aatlMd  1b 
tha  ^>aolal  oatagoiy  of  aaoond  ordar  dlffaraatial  aqaatioaa  la  vliiah  tha 
flrat  darlvatlva  doaa  not  i^ypaar.  A  ooa|ratar  prograa  la  daralopad 
uhloh  ganarataa  a  ntnbar  of  raouralon  f  onmlaa,  appllaa  than  aaaaaaalTa- 
ly  to  aolra  tha  initial  Taliia  problaa  7*H  7*0,  7O  ■  1,  7*(0)  ■  0, 
and  produoaa  sriqpha  of  aolutlon  arror  Ta  msbar  of  Itaratleaa.  Hm  ha- 
de prograa  la  than  aodlflad  to  aolra  tha  Couloah  Wcra  Ituatlea  •«''-*• 
<^(  /.  izr -  !di±})  : o    vlth  initial  oondlUona  app(rop(rlata  to  tha  aalaatad 
valuaa  of  tha  paraaatara  y  and  L.  Iha  Runga-Iutta  aathod  ia  uaad  to 
ooiqmta  al^t  valuaa  of  tha  Indapandant  Tarlabla  to  atart  tha  aolvtiaa* 
FlaaUjr,  tha  raaolta  of  thaaa  ooaputatlona  ara  dlaouaaad  In  datall, 

Soaa  of  tha  advantagaa  of  thla  aathod  orar  othara  la  eoaaon  uaa  arat 

a.  A  rational  oonvargant  la  ganarally  a  battar  apprcodaatlon 
to  a  fmotlon  than  ia  tha  oorraqpondlag  trunoatad  aarlaa.  Thla  vUl  ba 
dlaouaaad  furthar  In  tha  following  aaotion. 

b.  Only  about  ana  half  tha  total  nwbar  of  tabular  dlff aav 
anoaa  raqulrad  tagr  aathoda  uaing  polynoalal  appraadaatiflaa  ar*  naadad 
in  this  aathod  for  ooaparabla  aoeuraoy.  Ihla  la  baoauaa  dlffaraaoaa 
ara  dlatrlbutad  orar  both  aldas  of  tha  raouralon  f  onmla. 

o.  Sinoa  dlffaranolng  la  divldad  batraan  both  aidaa  of  tha 


•quAtion,  MitlMT  aid*  1«  diffiri«>i  too  mh^  Umo.     Ihia  n^ni  tt* 
nuiriE>or  of  ooRpaiotlflns  roqolrad  ad,  oon— q— itl/t  tiM  rMnd  off  orvor. 
d.    RoewdoB  fflCvolAo  darivod  tr^m  fmtm  ■iwfWfita  oro  •■- 
trwMly  o*i7  to  cMMMto  In  •  dlfiUl  M^«l»  fgugna.     ma  «1U  W 
lUustratod  in  SMtlon  <». 

2.      DrlTfttloD  of  Rocunloo  Fonwlai. 

AippoM  It  it  dosirod  to  oolTo  tho  initiol  ralo*  problaB 

y(x)  -  f(x.y(x)),      y(x.)-y.,      /•(«.) -y;.  (1) 

IntogTAtlAg  f  onuCLly  on  oooh  sldo 


•i.   A  '•   '* 

got  en  tho  laft  «ido 


Zhtroducing  a  cnango  of  VArloblo 

^  ■  x^  ♦  h»j^,      d<^  «  hd«^ 

t  ■  x^  *»•  hs, ,      dt  «  hds, 
and  defining  tho  forward  shift  oporator  E  auoh  that  K'"f(x)  ■  f(x««h), 
wo  get  on  tho  right  (treating  S  aa  an  algebraic  ^m^boDT 

'"2  ^2 


finaUy, 


.^,)V''-^'"'1^' 


or. 


2 


!  l^»  ' 


In  terms  of  th«  baokward  difr«r«QO«  opwator  v,  d«f  la«d  ty  v  /,*  ■ 
7;   -  7.-.,  .  this  is 

Multlpljlng  throui^  bgr  v*' 


^^Ing  throui^  bgr  V  i. 

slno*  V  operating  on  a  oonatant  ylalda  saro.  ^      ^ 

For  slapLioltgr  of  notation  «•  lot  s^  ■  ▼  7t  and  writ* 

Lat  a.'f  a,T  •«-  a^T  -f  ...  bo  tho  powwr  aarioa  aoqpanaioa  of  S  ~ — (• 

Tho  poNor  sorias  oan  bo  approodaatad  bj  a  rational  f^aotioB  bgr  um  «f 

tho  folloMlng  teofanlqxio.     a^-f  a,T  *  a^rS-  •••  <f  a^T*2    Ltfciftltif '<•"#»■ 

Multiplying  oaoh  aido  bgr  tho  danoainator  of  tho  abovo  f  roetion  and 

oquating  oooffioionts  of  tho  liko  powor  tonu  in  t  «•  ganM>ata  p*^i. 

liaoar  aquations  in  tho  f*q^l  wiknosna  b. ...  l^  ,  o  ...o^. 

b.-  a. 

b,  -  a,  ♦  a,  o. 

b,-  a,*  a.o,*  a.0^ 

b^»  a,-f  a,o.4.  a,  0,4.  a,o^ 

Iho  solution  of  this  qrstan  giTos  us  tho  eooffieiants  of  tho  r«t 

al  fraotion  idiioh  agrooc  with  tho  sorios  out  to  Fiq^l  tsms. 

It  should  bo  notod  that  a  giTon  polynoadLal  of  dogroo  r  eaa  bo  oquatod 


-A-V      I      V«* «»  ♦.  ' 


ld«ntlo«llj  to  Moh  of  r  ♦  1  di«^lIlOl  rrM^xons. 

n<»lal  «ft«r  1,  2,  ...,  n  t«ns»  MM«««iT«ljr,  w  obtAia  4  oollMU«n 
of  poijnot^Lal  fraoUons  irtilob  oaa  bo  aiTMicod  m  a  Pad*  Ublo  of  MtlMl- 
al  conTorg«t«.     Fig.  I  is  tho  Pado  tablo  asaoeiat«l  with  l^^']  i 
irhoro  a  third  dogroo  poljnoaial  la  naod  to  ganorato  tha  ratieaal 
oonrarganta*  . 

Iha  ponar  aariaa  «q>aaalon  of  )  — ^ \      ,  tr«neat«d  to 


torma,  la  tho  alxth  dagrao  poljooalal 


Iha  ooafflolMita  of  thla  poljnoalal  iiara  uaad  aa  lapvta  ta 
B  XSD  PADS  (Hayal  Poatgradoata  Sohool  Llbraxy  Routlna)  to  obtala 
coaff lolanta  of  tha  26  Pada  eonrarganta  nhoaa  niaMratora  «m1  daMitaatara 
ara  poljnoHlala  of  dagraa  q  and  p,  raapaotlTalj,  and  vhara  ptq  5  6. 
Pig.  2  llata  thaaa  ooaff lolflcta,  rotudad  to  thraa  algitfloatit  flforaa. 

Aa  atatad  abo^ra,  aaoh  of  thaaa  oonrargaata  ia  IdantiaaUy  aqtMl  to 
(3)  out  to  pfq4>l  tarsa.     Horaorar,  a  glran  oonrargant  opandad  algabrala- 
ally  Into  a  powar  aariaa  gaoaraUy  oontlnxMa  to  agraa  qvita  vail  with 
tha  ralatad  Inflnlta  avias  for  a  atabar  of  tanw  bagroad  tha  (ptq^i)    t 
l.a.,  a  rational  oonvargant  ganarally  ia  a  battar  approodaatlon  to  a 
function  than  la  tha  trunoatad  aariaa  froi  vhioh  tha  oonvargant  waa 
darlvad. 

Daalffiata  a  particular  oonvargant  to  (       ^        )      ^ 

Siq)poaa  wa  wish  to  find  an  ^>proaclBata  nuMrloal  aolutlon  of  tha  aquation 
y"  »  ^(x»y).     Subatltuting  (4)  for  J  -^ —  J    ia  aquation  (2)  «a  hara 

Rawrlttan,  thla  bacotaaa      {t^tJit^*   *ttyty,'iv')ns\^(f^^  »*^r%  . 


Noting  th*t  ^'U--  %:'^'''^    .  y^  •«  «P^"  «»•  t-^  •l^*i« 

Dmot«  the  oorffioi«it«  la  (5)  *•  foUoirsi 

Siibstltutlng  the  co«fflol«nt«  In  (5)  *•  hwret 
Q(0)«;-  QdK..*  ...J  Q(q)«,^-  h*[p(0)fi-P(lK..*  ...i  P(p)tJ)  .  or, 

«;-^Q(i)«i-.  -  ...  *  Q(q)«M  -^  h^fp(o)f^  -  ...f  P(p)f,-.]5  .  (^) 

Here  ei  is  given  in  terras  of  ooeffloienta  of  the  P«d4  ooPTergent 
and  of  previously  ooaputed  values  of  «  «nd  f .     Recalling  that  •;•  ▼  y,» 
Ji  -  ^,..  ■•■  y..»  we  miy  now  ealcuUte  y<  «  «;♦  ^^^  -  7,.,  •     **»•  ▼•^^  ^ 
f  •   U8«l  iB  (6)  win  h«Te  been  oon)uted  fron  a  first  app«al«ation  to 
7'  generated  by  a  predietor  fomula. 

3.       stability  Considerations. 

The  recursion  fomula  generated  from  a  given  Psde  ooRvergmt  is  of 

the  font:     Q(0)f,«  Q(i)»,. 1  Q(q)*.^*  h^[p(0)f;  -  ...  -^P(p)t  ,]>  (7) 

where  Q(i)  and  P(l)  are  as  defined  in  Section  2. 
Substituting  vV.«  y/-  2>v-*  yVv  ^<»i.in  (7)  »•  obtain 

Q(o)  (y,-^,,  *yis  )  -Q(0  (y«^ -^,.»  •V<.») -^^'i^  ^♦-f 

-^i^.,-*7iy.,)  ^  h^P(0)f^-  P(l)f<-.*  ...t  P(p)fw']. 

Multiplying  and  ooUeetlng  terms,  this  becoaes 

Q(0)y,   -  [q(1)  •»•  2Q(0)]  y,^   -  [q(2)  ♦  25(1)  ♦  Q(0)]  y^.^*... 

♦  (-1)*"  [Q(q)  *  2Q(q-i)  *  QC<i-2)]  yx-y  ♦  C-^)^  [«<<!>  ♦  ^<*-^>]  '^-r - 
+  (.1)^"  Q(q)y.^.,^**[p(o)A  -  ...♦  (-l/Kp)f,^"]  <•> 

The  general  eoluUon  of  thl«  difference  equaUon  Is  of  the  fe«« 


y^«  e.  i^-f  c,r^4  ...  +  o^p^i<i«r»  r, ,r«^««.,  r^ 

AT*  root*  of  th»  assocUtad  oh«r«ot«rlstio  •quatloo  (to  b«   i»tr»d»e«d 

lAtor)  and  tho  Cj«ro  oonst«nta« 

Con«ldop  tho  2nd  ordor  P«do  ooaTorgcnt       ,_^  ^  ^  ^^       (F«n»iU  13 
in  Pig.  ?),     Pron  tho  dlooussian  of  Sootlon  2  wo  knOTr  that  tho  roe«r«i«» 
formilo  ganorotod  by  this  oonrorgont  will  haTo  tho  followlnf  ooafTlolantst 
Q(0)  •  1;     Q(l)  «  0,  i  *  1 

P(0)  -4;  P(l)  -  -tf  I  P(2)  -5»  P(l)  -  0.  i  *  3 
tho  gonoral  roourslon  foranilo  (8)  than  boooaoa 

y,  -  ^... -  y..^*  ^t,  ♦/Jfw.  ♦4f...]  .  (9) 

Jn  tho  Initial  valuo  problaa  (1)  it  i*  ooaa  that  fi «  -Jj  • 
Substituting  in  (9)  and  oollooting  tormsi 

OP, 

(I2fh*)y,-(2iul0li*)y,..*  (12»h'  )y,^«  0 

Plndlng  root.,    /v  =   (x..-,on  ±  y  ( l4^u<.>>'^4.f  .a.<>;Sl 

(10) 


Ignoring  torms  of  tho  ordor  of  h  wo  haro 

'a. 
Tho  ganoral  solution  of  (9)  Is  of  tho  font 


(U) 


y^oro.  !;♦  o,r» 

Wo  dotonOno  tho  oooffioiants  o;  fron  th*  ■olutlon  y  •  ooa  *  of  th* 

initial  raluo  problon. 

y^  ■  e,  r  •f  o^r*«  1  (Initial  ooaditloaa  glv*  y.  ■  1) 

y^  •  o.  r>  o,^'"  1  -—      (y  »  cos  h  ■  l-i*!,,,,  ,  For  —11 

^        valuoa  of  h,i-^*ls  a  good  approila»tlon) 


Solring  thm—  •quationii  slaoltaiMoaalji 


Thus  w*  h«v*t 

V -- i^iiru**^')" -^  ii^p-».>r         (12) 

This  oan  b«  irritUn  In  trisonowtrio  fonit 


(13) 


for  auU  h,  h?#^lii  vhioh  oam 

and  this  aolution  ooafXfa  to  th*  aacact  aoliitloB. 

It  oan  b«  80«n  fros  th«  graphs  that  Fonmlaa  10,  11,  12,  15,  22, 
and  25  ara  tha  moat  proalalng  froa  tha  standpoint  of  alninRai  arrer 
•oplifloatlon.  Qnfortunataly,  thalr  Mssoolatad  oharaotarlatie  aqvatlaia 
(fourth  through  savsnth  dagraa)  ara  too  Intraotabla  to  paradt  analysis. 
It  would  ba  quita  siopla,  hovavar,  to  obtain  tha  roota  of  thaaa 
aquatlona  with  an  appropriate  oo^putar  prograa.  Tha  arror  oharaotaristios 
of  tha  aquations  oould  thsn  ba  datandnad  froa  a  stu^jr  of  tha  propartiaa 
of  thalr  roots*   "nd.!  is  dlsoussad  in  Saetlon  7  as  a  prolaot  i  anit— inlil 
for  futura  study. 

Tha  foragoing  pu-agrapni  wara  ooneamad  vith  tha  aootanua^on  of 
trunoation  arror  without  ragard  for  tha  faot  that  arrors  in  j^.,  and  j^.j. 
in  (9)  will  affaot  srdbsaqusnt  arrora.  To  oonaidar  this  sffaat  a 
dlf  f  ar«nt  approaoh  will  ba  taken. 


In  tama  of  th«  differanc*  •qoatlon  (9)» 

where  7,  is  the  ▼alue  of  solution  of  the  Initial  ralue  problea  at  x^ 
and  f^   =  f(x,  y, ).     T(h)  ie  the  tnawatlon  error  iBtradhwd  \^  wdsf 
,~^  to  aiiproxlJMte  f-j/'"^M     .  "id  R(i)  la  the  roondoff  mn». 
We  define  I(h,i)  «  T(h)  ■«■  R(i). 

Defining  e,  «  y;  -  y,  ,  we  subtract  (9)  fro«  (14)  and  obtain 

Asmaalng  that  f^(xi,y)i»  dlff»«rentlable  ulth  reapeot  to  y  la 
the  interval  (y^  ,y^  )  we  nay  uae  the  mean  Talue  theorea  to  writes 

^  ^^(^Ifii)  do««  not  vaxy  greatly  we  matj  let  fj  (x„<,.)  •  K  and  write 

Collecting  tenut    Cz-T^KV^t,-  f  Af  J^KV^tf^..  -  Ci-^A^*^*,.,  ♦fCrO 
or,      03-KV)<,  rra*4.ipkV^6^^   .  (la-iV^^.^    ^-  f  ri.«>V 

to  the  aolution  of  the  initial  value  problM  (1),  K  ■  f^(xi,ii)  •  •!• 
Hence,     (12+h'^  )«4-  (2iul0h*  )e,^  -  (i2«i*  )«^.+  E(h,l).  (15) 

Itoder  the  a8suiif)tian  that  S(h,l)  li  r*laUvely  iaaanaiUTe  with 

reapeot  to  1,  we  niy  proceed  as  ve  did  in  aolTlng  (9).     (Aa  implied  to 

this  oase  the  truncaUoo  error  T(h)  «  Ch*y  *^(i>   ,  and  j'**(t,)  Tariea 

slowly  whan  the  h  incr«Mnts  are  muH).     We  will  h«nceforth  ai 

E(h,i)  to  be  a  constants     Eh^. 

We  Host  sdTe      (12  ♦  h*)#i  «     24-lOh'  )«,,  -  (I2fh*  )*^  ,♦  Dh* 
«•  «  0  (no  error  inltialiy),  asd 

«•    ■  /*         (aasuMd  to  be  non-^ere) 


Ibe  soldtion  oonsists  of  the  general  solutioa  nnnTSipwlli^  to  th« 
hoMganeous  part  ^«  0,  r7*  e«r/  and  a  particular  aelvfeloc 


ponding  to  th«  non-haaoc«MOoa  p«rx  un  .     tn«  p«Lru.oiu«r  wxaiion  9mk  xtm 
found  bgr  asmadng  thatC^la  a  eonatant  A,  ao  that 

(12fh*  )A  «  (2^10h'  )A  -  (12^^  )A  -*-  Ch^ .     BaaM, 

A  r  PK^ _     _oi^t    -      2iL  • 

Tha  o(aq>lata  aolutioa  la  than^«  e,  rf<f  ^r»  -f^^h  .  (16) 

Lat  r,  ■  a  4  lb,  r^«  a  -  lb. 

Dmota  o,  «A4-  1^,  o^-Jl  -  1^,  whara  -f«  -  1  h^'and  ^: -^  (/♦-(•  j>iK*|  * 

Tha  aolution  (16)  la  tranaforaad  Into  trlconaaatrle  fom  b^  lattlag 
&  >  arc  Bin  * 

Than,e^<B   |rl   (2Aooa  nO  -  2^aln  n6)  ♦  Jit  h",  nqpaadlagt 
*.-  |rr^-^h*^ooa  ne  +^(/-»'  (a-l);|h*)  aln  ne    ♦^b\ 
In  partloular,  if  wa  approximta 

a  •  1    b  ■  h  ^  ^'^•f^M   al-t-lh,  r^al-lh 
6»  »  (fT+K"*^ )"/  -  j^  h'  ooa  no  ♦^C  /aln  n«)'V  -i-;!  h^  (17) 

wfaara  9  »  aro  aln  ^^-1»     .     For  aaall  h 

(i8) 


'i*k 


6^  i-  (ttf?  )"  >  -  h*'oc8  nh  -f/ain  nh  V   ■••—  h*. 

Thla  aolution  of  the  arror  propafatlon  aquation  quail t«tiTalj 
plains  tha  baha^rlor  of  tha  flaq>arliMntal  raaulta  diqtiajad  in  fif .  16, 

a.  If  h  la  aMdl,  tha  alna  tam  of  I  .±  h^eoa  nh  4^Kln  nhl 
la  aignif leant  and  tha  eoalna  tam  la  naglifibla.     brer  would  ba 
axpaotad  to  Taxy  naarlj  ainuaoldallj. 

b.  6.^  ihould  ba  approximataly  aqual  to  taro  lAaira  nh  la  anl- 
tipla  of  Tf,  dnoa  tha  tana  -^  aln  nh  vlU  Tanlah  at  thaaa  polnta.  (UTan 
h  >  .02  and  aolTlng  for  n  «   J^    ,  k  ■  1,  31,  ••.,  5»  va  aoqpaet  to  find 


nmr~m*TQ  •rror  «h«r«  n  ■  157*  3i^»  ^t  6>B  iad  785*  Flc*  16  ^mn  tlila 
to  b«  th«  OAM.  Error  nulla  riioMii  ia  fie*.  33  vd  ^1  ••■  b* 

■lailarlj  bgr  aolvlng  for  n  with  h  «  .Ol  aid  .05t  raapaotlTaly. 

o.  Th*  braokotod  torn  In  (18)  ia  ■altlpUod  b^lrl"!  imd^, 
tha  arror  Bagnltuda  ahould  grow  with  ineraaalac  n  If  ill  >  1  and  dlmiaiah 
IflrUl.  Tha  anplificatlon  factor  {.yn^*f   in  (16)  ia  onlj  m 
■ation  of  Irl.  Iha  actual  roota  i*ara  foond  to  ba 

for  all  h,  Uhanavar  h  ia  aaall  anon^  to  Inaura  that  (16)  aloaalj 
laataa  (17)t  arror  naxiaa  ahould  raaain  naarlj  ccnatant.  Fie*  33  a«b- 
stantlatas  thia  for  h  »  .01.  Ihcraaaing  arror  aasdjaa  in  Figa.  19  nd  kl 
(h  ■  .02  and  .05(  raapactiraly)  ara  attributad  to  tha  failura  of  h  to  aaat 
tha  oritarlon  atatad  abova. 

No  axplanation  ia  offarad  for  tha  hi^  fraquao^  arror  oaeillatioa 
ao  pranralant  in  aoaa  of  tha  grapha,  a.g.  fig.  16.  Unraatigation  of  thia 
phanoMnon  is  suggaatad  aa  a  projaot  for  futura  atud^. 

4.   Paeorlptlon  of  tha  Taat  Prograa. 

Thia  prograa  ia  daaignad  to  aolva  a  giT«i  aaoond  ordar  diffaran- 
tial  aquation  26  tiaaa,  auooaaaiTaljr  uaing  aach  of  26  Pad^  aomrarganta 
to  ganarata  a  oontral  racuraion  fonnila.     Litagration  of  tha  giT«i 
aquation  ia  by  tha  pradlctor-oorractor  aathod,  ualng  tha  Mixta  pradietor 
foraulai    y,  y..  *.^,.,  -  y  .^  *!*  jV,:.,  vi,:^  f^^.^]  -  laaaaaaiy 

atarting  valuaa  ara  ooiqf>utad  bgr  tha  Runga^utta  aathod.     Iha  itairatiTa 
process  may  ba  oarrad  out  to  aa  aaxgr  aa  900  stapa,   aolution  arror  baiag 
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d«t«nitn«l  At  MMh  at^  If  «n  ouwt  •alwtiaft  of  tb*  giT«i  «qp*tiMi  la 
known.     Qr^ihs  ar«  plottad  with  •»«  *■  Mdiaat*  wd  aaribar  af  11 
«a  abaolaaa.     rig.  3  axplAlna  tha  ^iibela  tsMad  la  tha  prcgna*.    Fig,  ^ 
la  a  printout  of  tha  ptognm  aaad  to  aralvata  tha  aquatian  gaaavAtlag 
tha  ooalna  function,  7"  4-  7  ■  0  (Progri  MOD  2).     Flf.  5  !•  *  prlBtavtt 
of  tha  procraa  uaad  to  araluata  tha  Couloab  Wa:Ta  Itiuatlan   J'4^(l'^'*^^tei, 
(Program  TMOD  6)* 

5.   Dlacuaalon  of  Data. 

a.  Prograa  IMCO  2 

(1)  Plga.  7  throng  25  •>«  cn^ha  of  pgogr— alTa  arror  ra— IW 
Ing  whan  tha  Initial  ralua  prohlaa  j*  -f  7  -  0,  7(0)  •  1,  7*(0)  ■  0  la 
aolrad  using  an  Intagratlon  atap  langth  of  .02  and  eorraetor  fonnlas 
darlrad  from  tha  rarloua  ?$di  oonrarganta.  Pig.  6  ralataa  tha  raauraioa 
fonnila  nuabar  Aam  on  aaeh  gnqph  to  a  daaorlptlon  of  tha  eorraapoodlBg 
Pada  fraotlon.  Porvolaa  16,  17*  20,  21,  23,  2i>  and  26  dlrairgad  rapidOj 
and  no  aaaoolatad  grapha  vara  plottad. 

(a)  Ponrolaa  10,  11,  12,  15,  22  and  25  produead  r«MTkabla 
aooura^  orar  tha  ontlra  raaga  of  900  Itaratloaa. 

(b)  Porwila  6  (Fig.  12),  although  laaa  aeeurata  th«i  tha 
gro\9  Just  dlao\isaad,  haa  an  lsq;>ortant  propartji  arror  aaadiM  ara 
dlalnlahad  vlth  Ineraaalng  Itaratlon. 

(o)  Tha  ranalnlng  grapha  Indloata  that  arroir  aaxlHa  ara 
not  balng  daapad  but  ara  Inoraaslng  llnaarly  vlth  raqpaet  to  n,  tha 
numbar  of  Itaratlona. 

(2)  Tha  naxt  aat  of  data  vaa  obtalnad  b7  rapaatlng  tha  pr»- 
oadlng  aoqiarljMnt  aftar  daoraaslng  tha  Intagratlon  st«p  langth  to  .01. 
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Th»  following  w«r«  obs»rT«di 

(a)  Thm   runs  uhloh  prorad  to  b«  ao^  aoetirat*  with  •  aU/p 
length  of  .02  (Fomul**  10,  11,  12,  15»  22  and  25)  w»r«  found  to  b* 
lose  aoourato  with  tho  aullar  stop  lanfth  (eoDtr«i7  to  what  would  b« 
«xp«jt«i),  FomuU*  15»  22  and  25  ahow  a  Mkll  orror  a^pUfieaUoo  iriiila 
all  tha  othara  In  this  group  Aow  a  aaoh  gr—Uir  arror  a^ilifieatien. 
rigs.  28  thrci;gh  33  ara  graphs  ihowing  tbasa  rasults. 

(b)  Femila  6  prorldas  1*m  aoourata  raaolts  than  with  th« 
graatar  stop  l«igth  bat  its  error  damping  oharaotazlstias  raaaia.  (rig.  27) 

(c)  Fonmla  5  is  aora  stabla.  Its  aoeuraqj  Is  oa^arabla 
to  that  in  tha  praoadlng  run.  (Fig,  26) 

(d)  All  othar  fomulas  prorad  to  ba  aora  aoourata  Uun  whan 
uaad  with  tha  graatar  stop  l«ngth.  All  of  than,  howvrar,  amplify  arrors. 

(a)  Tha  arror  graphs  assooiatad  with  Fomolas  9>  10,  11, 
13.  1^.  18,  22  and  25  vara  idantioal.  (Fig.  33) 

(3)  Incraaslng  tha  stop  Isngth  to  .05  had  tha  following  affaotsi 
(a)  Tha  fomulas  wfaioh  prorad  bast  with  a  step  langth  of 
.02  ara  much  laaa  aoourata  with  a  stap  langth  of  .05.  Figs.  3%  through 
M  ara  tha  arror  graphs  assooiatad  with  Formilas  6,  10,  11,  12,  13»  15* 
22  and  25  using  tha  .05  stap  langth.  Thaaa  gr^s  show  tha  on]j  foraulas 
rataining  any  arror  damping  oharaotarirtios  to  ba  nuabars  11,  12  and  25. 
b.  Progran  TMQD  6 

(1)  Tha  CoxJomb  Wcva  Equation  was  solrad  for  two  eonbinatlans 
of  Taluas  of  L  •BdY(l^l»7*i{  »*>d  1^2, 'y«3)*  and  with  various  inta- 
gration  stap  langths.  It  ma  not  f aasibLa  to  ganarata  arrer  |,mlia  far 
thasa  aquations.  Fig.  42  Hats  nagnituda  of  dariatlon  of  eo^otad  aal* 
utiona  from  trua  valuas  aftar  vha  nuabar  of  itarations  indieatad.  Dua 
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to  Ujnltatians  of  the  t«bl«8  uMdm,  th«  alrUnH  •rror  dlaovriMtal*  wm 
IE.7  (1  X  10'    ).     An  wcanination  of  Pl«.  >*Z  rmr—ls  tha  foUowlmt 

(a)  Id«rttloal  r«mlta  tmf  obtaliMd  in  all  •«•••  liMr* 
FormulAs  9,  13,  l'^,  15  >nd  18  ««r«  uMd.     This  groi^  also  produoad  tlia 
■oat  aeourata  aolutlona.     Radioing  atap  l«icth  laprcnrad  aaouraoy  throvgli- 
out  tha  ranga. 

(b)  Ponmila  8  aoluliona  -mta^  allchtlj  l«aa  aoourata  than 
thosa  of  tha  group  abonra  In  solving  tha  aquation  7"  -f  7  (1*-*1)  ■  0. 
Solutions  of  tha  aquation  7"  -f  7(1--  -■^)  *  0,  uaing  a  stap  laogth 

of  .02  vara  aora  acoxirata  than  with  anj  othar  fonralt. 

6,     Conclugjons. 

In  his  NiBarical  Mathod8[33>  Kopal  illuatratad  tha  ■hortooBlngs  of 
golng-ahaad  raoursion  fonnilast     thosa  darivad  froB  tha  ralationahip 

tainad  bafora  7^-  oan  ba  oaloulatad,     Tharafora,  tha  aathod  suggaatad  tqr 
Kopal  is  asaantiaUy  tha  pradlotor-eorraotor  mathod  ooi^tlad  with  a  pro- 
oadora  for  darivlng  an  unliaitad  aappij  of  eorraotor  f onnilaa«     9aaa  of 
thasa  wa  hara  saan  produoa  quita  aoovrata  aolutions  to  sararal  dlffarvA- 
tial  aquations  of  saoond  ordar.     A  glanoa  at  Figs.    '^3  -  ^5»  tha  mrrvr 
of  7"  4-  7  ■  0  aftar  aolutlon  by  tha  MUna  aathod,  la  aridooea  of  tha 
iority  of  a  nunbar  of  tha  oorraotors  darivad  froa  Pada  oonrarganta. 
It  is  suggastad  at  this  point  that  tha  lass  aoourata  goiag-ahaad 


typa  formulas  night  ba  gonaratad  in  tha  aaaa  wigr  and  uaad  aa  pradietors* 
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d«riT«d  f ran  a  pdljiiaHial  appraadaatian  to K  — ^^  \      af    imm   wi^i  n^ 


Ita*  faet  that  oartaln  oorraotora  dwrlTad  tram  Pad4 

tbair  ta#ca  ao  wall  riumLd  prorlda  aoff laiant  aotlTatiaA  far  a 

•tw^  of  pradjotor  darlTation. 

lona  of  tha  oorraotor  fonmlaa  diaevaaad  in  thia  papar  ara  atafcla 
for  all  aiiuatlcBa  anr  all  rang**  ft  ^     iddltleaal  fonaalaa  aaaM  ba 

but  thla  araa  doaa  not  appaar  proalainf  .     In  daallnf  with  aaaeaiatad 
oharaotarlatic  aquations  of  hi^ar  dmgr—  it  will  ba  a»tr<— 1/  diffla«lt 
to  looata  Taluas  of  paraaatara  whloh  will  ratala  all  roota  witlda  tiM 
raquirad  bounds. 

Tha  anallant  raaulta  obtainad  with  TotwkIm.  15  (Fif*  21)  iaAlMto 
that  a  fourth  dagraa  polomoid.al  approodMUon  to  i  Zlliil!?  V  la  aafflalait 
to  inaura  that  truneation  arror  ia  not  a  aignfieant  prebl«b     It  woold 
wmkf  tharafora,  tliat  tha  uaa  of  poljnawiala  of  hifhar  dagraa  waald  la- 
▼ita  additional  atability  eo^plloationa. 

In  auHBaxy,  it  ia  oonoludad  that  tha  proaadura  auggaatad  tgr  Kopal 
riiould  ba  Tiawad  aa  a  nathod  of  darirlng  pradiotor  and  oorraotar  f  anailaa. 
Tha  raaulta  aohiarad  with  aoaa  of  tha  oorraetora  darirad  in  thia  naMMr 
indieata  that  thia  prooadura  ia  praaialng  and  narlta  furthar  atvi^. 

?•       Bai?  niwamiations  for  Futura  Studiaa. 

a.     Invaatigata  tha  roota  of  tha  oharaotariatio  aquatifloa  aaaaaiatad 
with  tha  nora  prowl  aJng  raoardoB  fonmlaa  dlseuaaad  la  thia  papar.    A^ 
oartala  and  tabnlata  Talaaa  at  (  far  whiah  all  roota  Ua  wtthia  tha  unit 
olrela.    Vflth  thia  data  arailabla,  oaa  oould  raadily  aalaat  tha  baat 
raei'raion  f oivula  for  aolTing  a  giTan  diff arantial  aqiMktiea  orar  a 


ik 


sDoelTisd  rane*  ot  Taluas  of  th*  iadaponknt  Tarlabl*. 

b.     Dsrlve  a  Pad*  t«bl«  of  rmtlOMl  owwmiU  U^  (•••)'«^ 
f  ron  whloh  to  obtain  a  nwib«r  of  golag-alMad  r*eurslen  fomilaa.     9W>- 
■tltuto  oaoh  of  th«s«  In  turn  for  tho  pradlotor  f omnia  «a«d  in  Fr»fri 
T.-IOD  6.     Caloulate  errors  and  oooparo  with  Pig.  t^Z,     In  thia  ■iiar , 
evaluata  the  marlta  of  prwliotor  f oraulaa  darlTod  from  rational  Pate 
eonvargmts.     Ralston  ^1  dasoribas  the  adrantagaa  of  utlac  pradiotora 
and  corraotors  of  tha  aama  ordar. 

e.     Bictand  tha  poljmomial  raprasantation  of  |  — ■ T     b^rand 

ser^an  tarns  and  dariva  highar  ordar  raouralon  fonnalaa.     iBVMitlgata 
tha  roots  of  thasa  fomulas  as  Mggastad  Jn  subparagraph  a  abofra. 
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sD»cifl«d  rang*  of  TaluAS  of  tb*  liiiWp«nd«tt  varlAbl*. 

b.     DvriT*  a  Pad*  tabl*  of  ntlonal  oanrtr^tmtm  to)  (>-*>' 
f  ron  whleh  to  obtain  a  nunbor  of  golac-idioad  roounian  foHMlaa. 
stltute  aaoh  of  thoso  in  turn  for  tho  pradiotor  f oraol*  ttood  in 

Ti-KX)  6.     Calculate  arrors  and  ooi^aro  with  Tig,  42.     In  thia  MHMr« 

# 
waluata  the  aarlta  of  predictor  fomulas  dorired  froa  rational  Pad* 

convergants.     Ralston  [S]  doaoribes  the  adranta^ea  of  uaiiic  prodietoro 

and  correctors  of  the  saa*  order. 

0,  Ebctend  the  polomooial  repreaantation  of  i   — -    \   b^fond 

seven  terms  and  derive  higher  order  recursion  f  onnilas.  Jhrostigato 

the  roots  of  these  foinnvlas  aa  auggeated  in  subparafraph  a  above. 
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rig.  1  .  3«q>l«  Pad*  Tahl* 


V^ 

0 

1 

a. 

3 

0 

1 

I 

3 

1 

» 

'♦'♦V'^fr' 

i 

1 

1 

t 

«-^*.ir* 

'-'♦rir» 

1 

l-r^.;^r*' 

^^'f 


d«ilT«d  from  third 


C«RT«rs«nt«  to 
poljnaHial. 

Tbm  f  Mt  that  thr««  of  u  w  o«nr«rc«Bta 
•qual  la  oolneldaBtal. 
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Tlf.  2  -  Ce«drriol«nta  of  Pad*  ConwrnrfmU  to  i  "'^  ^'    '^  f 

Mo.  Iwmtor                                                D»o«ln«t«» 

1  1.0  1.0 

2  1.0.1.0  l.O 

3  1.0,1.0..917                                             1.0 
**  l.0,l,0,.917,.833                                  1.0 

5  1.0,1.0,.9l7,.833,.76l  i.o 

6  1.0,1.0,.917,.833,.761,.700  l.o 

7  ^-0                                                            1.0,-1.0 
^  ^•0-^3                                                  1.0,-.917 
9  1.0..09l,.008                                         1.0,-.909 

10  1.0,. 087,. 003,-.  00«*  1.0,-.913 

11  1.0,.080,-.003,-.010,-.005  1.0,-.920 

12  1.0,.07^-.009,..016,..011,-.005      l.0,-.9a6 

^^  ^•°  1.0,-1.0,. 083 

'^  ^•O'O-o  i.o,.i.o,.oe3 

15  1.0,.600..050  1.0,-.^,-.500 

16  1.0,-1.29,-.121,-.014  1.0,-2.29,1.25 

17  1.0,-.830,-.082,-.013,-.002  1.0,-1.83,. 832 

^®  ^'^  i.o,-i.o,.oe3 

"^  '-^-^^^  1.0,-.917,.003,.00«» 

20  1.0,-1.0.. 050  1.0,-2.0.1.13,-.133 

21  1.0,..95^..077..002  1.0,-1.93.1.12,-.135 

"  ^•°  i.o,-i.o,.oe3,o.o,-.oo<» 

^^  i-o.-i.o  i.o,-2.o,i.oe,-.oe3,-.oo* 

2^  1.0.-1.0,.l23  1.0.-2.0,1.21,..206,.006 

^^  ^'^  i.o..i.o,.oe3,o.o,-.oo^.-.oo* 

^  ^-0.-877  1.0,-1.88,.960,-.073,-.00k,-.001 
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rig.  3  -  L»g«nd  of  FTogrwi  4n*ols 

X  Ind«p«Kl«nt  Tartafal* 

I  D«p<nd«nt  ▼ariabl* 

IP  y» 

TDP  y* 

IT  y 

EIUIGR  7-7 

D(I)  ^ 

z(i)  ▼'y 

Ali(2)  Cccffioi«nt  of  2*^  noMrator  tmcn  in  Pad*  e«nTWi«nt 

ADC^)  Co«fflol«fit  of  k      dmoalnator  torn  In  Pad*  oowrvrgaBt 

AX,  AT. 

AXF,  AIDP    Initial  Taluos 

JP  Niaibor  of  first  oonrorgvit  to  bo  oraluatod  (l<jptf26) 

JQ  Huaibor  of  laat  oonvorgant  to  bo  ovaloatod  (l'JF'JQ'26) 

H  h  (intogration  ct^  lancth) 

N  Rumbor  of  starting  raluos  to  bo  ooapotod  fcor  Kwgo-Kutt*  — tlKJ 

M  Total  niBBbor  of  itoratlons  dosirod 

BI(I)  Danotos  numbor  of  torma  in  nuMrator  of  fonnila  no.  fZ). 
Usod  for  information  in  printout 


ID(I)       Samo  aa  I1I(I)»  aacoopt  rafara  to  daneainator  torma 

DBL(K,I,J)   Array  vhioh  storoa  Pada  oooffioianta.  K"l  danetoa 
nuaorator  ooaffioiant;  1*2,  danoadnator  eooffioiant. 
I  ia  indaK  no.  of  eonrorgant.  J  ia  tha  mater  of  a 
apooifio  torn.  I.g.,  DKL(2,7i'^)  daaotoa  tiio  ooaffioiant 
d  tho  foarth  torn  in  th*  danoadnator  ot  aowr argent  ae.7* 


00  Cooffioiant  of  a  ^  in  roeuraion  fonnila 

Q(J)  Ooaffioiant  of  a,.;  in  roouraion  fonaola 

PO  Cooffioiant  of  y*^-  in  roouraion  fonmla 

P(J)  Cooffioiant  of  y*j   in  roeuraian  fonwla 

L  Ihdaac  no.  of  eonrorgant  and  ita  aaaooiatod  rootaaian  fotaala 

18 


Fig.  U  -  Provjc^a  IWOO  2  (Prlnt«»t) 


r  ^  o  •  •  A  T      f  1 


M  1  '^ ) 
!  1.  1  /»?Il.lf"Fl'^.]"> 
.M'^    H=,    F6,?,     V<         N-  =>H         M«»     14//) 

r 

"  r  »  ^     '  i"-  « 

!T  =  1  .1'' 
-     !TI  'ff  IT)='»^-' 

IT1TL«^(P)  =8MF,  T'r-I 
ITITl.^{9»=8H'"i.  '="RROR 
TTITL«^(  1.'^)=«H    V«^     T. 

^  r  !  '-  I  r  M  T  - 

"  .^,<,  I,(  IN(  n.  !n(  n  ♦('^fL(<*I»J)  •  ' 

f^T  I  Ml  If" 


«,  -  '> 

v 
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flj.  >*  (cont'd) 

QKDr  (Tif  (  ]  )4-      ' 

X(  I  I«^X*H 

V(  I  »=f^Y*V  +  '^<T 

Yr,p(  J  T  ,  , 

RXaX(  . 
P.V»Y<  I  ) 

V  'MT  -  V  "\  D    (     t    > 


r  rovpijTF  Fi?^T  ANO  ^^ronn  r>\F.froFnr^<,  or  y 

no  91    i»2,N 
^i  r>(  I  »=Y(  I  )-Y<  i-i ) 

00  92  I»3tN 
92  Z(I)=0(I)-D( l-l) 

* 

L«JP 

C 

c 


T/-\  n  r 


rrMPltT«^  COFTFiri'^NTS  OF  RrCllC*;  T'">N  FO?M»JL*  FPOM 
Cr^FrirfTNTS  OF  P^nF  ''ONV'^R'^-'^MT 


AMt 1  )  =  ]  ,0 

AO(i)=i.n 

2'^^  CONTINUE 

no  207  J=2f6 
AN(J)=DEL{1.L.J) 
207  A0( J)«0EL(2fL»J) 

12  OObO.O 
203  Pn.0.0 

00  14  J=l,6 
202  Of^sQO  +  ANfJ) 
1^  oospn+An(J» 
00  15  Jrl t^ 
0(  J)='^'.0 
P(J)=0.0 
00  15  I- J. 5 

0(  J)=0(  J)*AN(  I  +  l  )  »L;riO  •',  (  1  tJ  ) 
15  P(J)=P(J)*Ar*(I  +  l  )»=»INOM(  I.J) 
» 

NN«N+1 

DO  I')  1»NN»M 
* 

r       COMPUT-  APPPOXFMAT*"  VALUF  OF  Y  WITH  MILN^  PRFOICTOP 
Y(n=Y(I-l)+Y(l-3)-Y(I-4)  +  (H«»:»/^.'^»»(^.'^"^^''''-n 
1      ♦7,»Y0P(  I-?)*5.»YnP(  F-'<)  ) 
YnP( I )«-Y( I > 
• 

C      USF  RFCURSION  rOI^MULA  PFRIVFO  FROM  PAOF  CONVE'^GENTS 
C      TO  COMPUTE  SErONO  DIFFERFNCFS  OF  Y 

0R«-0( 1)»Z< I-l) +0(7) •Zf 1-2 >-Of*»*Z< 1-3 » ♦0141 
I     «Z( I-4»-0(5)»Z<I-5) 
PB«(P'^«YDP(  I  )-P(n«Y/>P(  I-1)+P(?»»Y0P«  I-2»-P(  3» 

Zl I )>IPB-OR)/00 

20 


on  )«z(  n*t)i  i-n 

19  Yr>P<  I  )»-Y<  I  I 

• 

C  COMPUTP    TRUF    VALIJF*;    OF    Y    AND    raLTULATF    FRBOR 

210    00    20    I«1»M 
XKD-I 
AI-I 
YT( I )«rOSFIH»AI ) 

20  ERRORS  I  »«Y(n-YT<  I  ) 

LABEL-tBCO(LtOf^*IER) 

CALL  DRAW  ( M*X t tFRRORvOtOtLAnPL* IT ITLEvO. tO* vOtOtOtO* 
1     5.5»ltLAST) 
IP  (L-JO)  24,25.2^ 

24  L»L*1 

GO  TO  206 

25  PRINT  999 
END 

FUNCTION  8IN0M  (M,N) 
AM>M 

IF  (N-0)301. 302.302 
302  IF  (N-l)3O3.3O4,305 

305  IF  (M-N)301,3C3»3O6 

306  CONTINUE 
IC»N-1 
NUMsM 
N^NOM.N 

DO  444  I»l»< 

NUM«NUM»(M-I ) 
444  NONOMxNDNOM»(N-n 

XNUM=NUM 

XD«NDNOM 

BINO.M«XNIJM/X0 

GO  TO  309 
•i0  3    BINOM.1,0 

GO  TO  -^99 
304  BINOMsAM 

GO  TO  399 
301  PPINT  312 
312  FOH'-'AT  M  fiH  RINOM  INPUT  FAUtTt 

399  EN 
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rig.  ^  -  ^'^^^igtm  MOD  6  (PrlaU«t) 


» 


♦AJIO' 


.'♦l  1'  ) 

f  IP2E15.7.0P2F10.'^»IPE15.  ^  i 
100?    FOHi^AJ     (4H.">    H««F4.2t*H       N»»1^»4M      M*«I4»%H    ETA>tF4«l« 

1  -,i\    L«»F^.l//) 

IDOA    F  (  Il»I?i2IUlP'iFl''.10) 

''^  '^  '     ..  ,  (laH."!  FORMULA  NO.  *\2*f>H»       NUM  •,I2» 
9H   DFNOV  x,I2//) 

1007  FOPMAT  (15H''  JOR  COMPLETED) 

C      READ  INPUT  DATA 

READ  1000»  JF.  JO 

READ  1001tH,N»M»I'-itIT 

READ  1002*  AX»AY,AYP»ETA»ZL»CO 

PRINT  IC'Sf  H»N»MtETAfZL 
» 

C      CONVERT  INITIAL  VAL'JFS  FROM  NORMAL  IZ'D 
C      FORM  USED  IN  MOST  TA»?LES  INTO  fTANDARD  FORM 
C      DFSCRIPrO  IN  TWIS  OROGRA^  AS  V 

LL«=2L 

CL  =  CO 

DO  201  L«1.LL 
201  rL=(CL»SORTF'ZL»*2+ETA««2)  »/ (7L*I  2«»2L+H  ) 

AY=AY»CL»(AX«*(LL+l» ) 

A(LL-»-H"l. 

A(LI +2)«ETA/(ZL+1. ) 

JL  lL+3 

DO  202  L«JL.100 

XL  =  L 
70  2  A(L)=(2.*ETA#A(L-1 »-AIL-?) ) / ( ( VL+'L I»( VL-ZL-l • 11 

PHI<^«'".'^ 

JL=LL+1 


20  3  r-Hl  %  =  PHlb  +  XL»'v  I  L)*AX»*(L-LL-1  ) 
AYP:=CL»PHIS»(AX*»LL)  . 
♦ 

C      READ  PAOE  COEFFICIENTS 
DO  11  L«l»52 
11  READ  10<^<*»  <«I  «(  INCI  )  tlOI  I  l*IDFLIlC*I*J)»J«2t6M 
» 
'-       ^OMPUTP  N  "^TARTINr,  VALUF'  °^  o..ur.c  ..'iitt.  MCTunr. 

VbH»AYP 

X(1^)»AX 

Y(IS)=AY 

YP(IS)«AYP 

YDP(IS)— AY»(l.-(2.»ETA)/AX-CZL*<ZL*l.n/AX««2)  ^ 

lSl-IS+1 
IS2-IS*2 
DO  105  I-ISltN 


Fig.  5  (oont'd) 


100  Bx«x(i-n 

•v»v(i-n 

TO  lo-^ 

iM  riX«X(  I-n+H/2. 

=^Y«Y(  r-l  )*V/2.*RK(  n/4. 
'•,0  TO  103 

102  «x*x(i-n+H 

PY  =  Y(  t-1  )+V  +  f^K(  3) 

103  9YDP  =  -BY*(  l.-(2.»PTA)/')X-(ZL»(ZL*l.»»/flX»»2) 

104  «*K{  J)«(M»*2/2.)»BYr>P 
PKT«(RK(  !)♦"<  (2)  +  «iKL('<»  I/"*. 

RlCP«(Rti:(  1  )+?.»P<(2)  +  2«»BKn)+PKI*»  »/3. 

X(  1  )«X(  i-n+H 

Y(  I  )=Y(  I-1)*V  +  3<T 

V=V+BKP 

YP( I )=V/H 

105  YDP(  I  )«-Y(  I  )»(  l.-(2.»ETA)/X(  I  )-(ZL»(7L-»-U)  )/X(l)««2) 
* 

C      COMPUT.«^  FIRCT  ANH  S^^rONH  OIFFrpcNrPS  OP  V 
no  106  I=ISltN 
1)6  D( r )»Y( I)-Y(I-l) 
no  107  I»IS2.N 

107  ?'(  1  )«n(  I  )-0(  I-l  > 
» 

C      SELECT  PAOE  CONVFRr.ENT  TO  BF  USED 

L  =  JP 
S 

C      COMPUTE  KECURSION  FORMULA  COEFFICIENTS  FROM  PAOE 
C      COEFFICIENTS 

AN(l)'L.O 

A0(1 »«1.0 

108  CONTINUE 

no  100   J=2,6 
AN(JJ»DEL( lfL»J) 

109  An(J)«DELl2tLtJ) 
OO-O.O 

P0«0.0 

DO  110  J»lt6 

Qn»00+AN(J) 

110  PO«P0+AD(J) 
00  111  J«I»? 
Q(j)=o.n 
P(J>«0,0   ^ 
DO  111  I"J»5 

Q(J)«0( J)+AN( I*1)»BIN0M<|,J) 

111  P(J)»P(  JI*An(|>l»#(iINOM<I»JI 
« 

NN-N+1 

DO  112  I-NN.M 

X(I)»Xn-l>+H 
» 
C      OBTAIN  FIRST  APPRo'xi  MAT  ION  TO  V   BY  USING  MILNE  PREOICTOff 

23 


Fig,  5  (oont'd) 


YnP(l>a-yfit«(i«-f:-«cTA|/xil)-< 

* 

C  „,„    ,  ,  "c     -    OF    SFCONP    ...-rr-RPHrF    Of    y 

r  li^fN^'    PC  RAT="0    FROM    PAr>e    CONVP«»GPNTS 

p«=(Pr>*YnP(  I  )-p(i)»YOP(  i-n+p(7j»Yr)Pi  i-?)-Pii)«Y0P||*9) 

1     ^-^f '.  )»YnP(  1  -4) -P(  S  )«YDP(  I-«i)  >  »(H»»2) 
-06)/0^ 

« 

C      COMPUTE  Y  FROM  DIFFERENCES  JUST  CALCULATED 
0(  I  >=Z(  I  )+D( I-l ) 
Y(  I  )=0(  n+Y<  I-l  ) 

•  i«-Y(  n«(u-(z.»ErA)/xin-(ZL«(ZL*i.n/xin»«2) 

C      CONVERT  ALL  COMPUTED  VALUES  0*^  Y  INTO 

C       NORMALIZFD  FORM  FOR  '''^;<daR  I  «;0M  WITH  «,TAND'\PD  TA9LFS 

00  ?'>0  I«NN»M 
200  Y(  n=Y(  I  )/(CL«(X(  I  »»«<ZL+1.1  )  ) 

PRINT  10n5»L.( 1N(L>.ID(L> ) 

PRINT  10n6,(Y(  n,I«lS»M.IT) 

IF(L-JO)  113tll4tll^ 
113  L=L+1 

GO  TO  108 
lU  PRINT  1007 

FNO 

FUNCTION  BINOM  (M»N) 

AM»M 

IF(N-0)3'^7,301,'^')1 
^01     IF    (N-1 )305. 306,302 
^02     IF     (.M-N»'»n7,303«3O'^ 
303    <=N-1 

N(JV  =  M 

--H 

'     1  =  1. < 
iM-I  ) 

T;'i  A     ,  ,  .  inv<*  I  w_  I  \ 

XNUM  =  ' 
XD»NC>'l'i  , 
BlNO.'^rXNUM/XD 

r,o  TO  30fl 

PIN0M«1,0 

GO  TO  308 
106  BINOMsA*' 

'--OTO  308 
vv..  POINT  3000 
3000  FORMAT  (18H  BINOM  INPUT  FAULT) 

■■..YOP 


J» 


rig.  6  .  RalAtlon  of  FonulA  Jtwtmn  to  Pad*  CaHrwcvata 

forwul*     D»gr««  of  Coenr^rgent  fW/D)     Formula     Dbctm  of  0— rtgit  (MfD) 
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Fig.  7  -  Fornttla  I  Error 
aUp  langth  ■  .02 
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Fig*   12  -  Poraula  6  Bkror 
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Fif*  13  -  FonRil«  7  BTTor 
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FLg.  15  -  Formula  9  ■rror 
9Up  Icncth  •  .M 
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Flf*   16  -  FbnMl*  10 
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fig.  19  -  romul*  13  Error 
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flC.   ?7  -  Formula  6  frror 
St«p  l«mcth  *  .01 
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St«p  l«neth  ■  .01 


at  >- 


h'TT'f^^-W^ 


Fir.  31  -  ."ornul*  13  fiirror 
5t«p  Isncth  ■  ,01 


y.  > 


50 


fl»^     7?  .  Forwila  22 

St«p  l«f«th  «  .01 


11^ 

ir 


51 


'•r? 


J} 


Step  lencth    •  ,91 


mA        '' 


13 
l««tii 


.01 


^     ! 


LJ 


CO 

CO 

Ill 


I 


n'.'-^jlJWp-,'^,' 


5t«p  Iwn.rth  =  .05 


CO 


lis 

t  •  >- 


X.  >- 


i 


53 


r 


fm*-^'T^^^r  j^,mvw-Tmi^rn^,w  wm^tw^^ 


^t«p  length  •  .05 


LiJ 
CO 

l-H 

CO 


12 


fig,  36  -  Fomul*  11  Srror 
9fp  length  •  .05 
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?ift.  T7  -  ?or«ul*  12  Brror 
3t«p  l«neth  =  .05 
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or'.ula  ?2  Error 
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fig.  42  .  Solution  orrort     Couloab  Vtcwv  liiiiAtlana 
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